Effect of weak magnetic field on quantum spin correlation in an s-wave superconductor 
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We study the effect of a weak magnetic field on the state of the spins of a pair of electrons forming a Cooper 
pair in an s-wave superconductor. By perturbatively solving time-dependent Bogoliubov equations up to first 
order, we obtain the two-particle Green function, and whence the two-electron spin-space density matrix. It 
appears that up to first order approximation, the spin state of a Cooper pair retains the form of a Werner state. 
This state is then examined in the sense of quantum correlations, i.e., quantum discord. 
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In quantum manybody or strongly correlated systems, 
quantum correlations are responsible for a host of interesting 
physical phenomena. Specifically, quantum entanglement is 
known to be at the heart of quantum phase transitions and a 
major player in quantum information science and its appli- 
cations (30 . It is thus of fundamental and practical impor- 
tance to analyze quantum correlations generated in quantum 
manybody systems. As an example, recently it was shown 
|0] that in a non-interacting electron gas and at low tempera- 
tures, spin-spin entanglement behaves quadratically with tem- 
perature. Additionally, two electron spins forming Cooper 
pairs in an s-wave superconductor has been shown to be a 
Werner state, whose entanglement may survive within a co- 
herence length £o [5]. In quantum correlated systems such as 
extended Hubbard model, whose states are amenable to ex- 
act or approximate calculations, the role of entanglement(i.e., 
quantum correlations) can be identified more explicitly ||6t]. 

In this Brief Report, we consider an s-wave superconductor 
subjected to a weak magnetic field. Of our special interest is 
the state of a Cooper pair in such a system and how apply- 
ing a field can alter quantum correlations between the spins 
of the electrons in this pair. To this end, we first obtain the 
density matrix of the system by employing the Green func- 
tions and solving the related Bogoliubov equations perturba- 
tively. It is shown that this results in a Werner form for the 
spin-space density matrix. Additionally, we calculate quan- 
tum correlations — exactly, quantum discord [7] — of this state, 
from whence one can see how applying a weak magnetic field 
modifies quantum correlations in the level of a Cooper pair. 

The Hamiltonian of a generic superconductor subjected to 
weak magnetic field is as follows BS MlOTI : 

H eS = J2 [ *tW [#e(r) + U(r)] ^(r) d 3 r (1) 

+ y"[A(r)*t(r)*j(r)+A*(r)^(r)* t (r)] d 3 r, 

where ^ (\&) is the fermionic field creation (annihilation) 
operator, H e (r) = ^[jfiV - eA(r)/c] 2 - fj,, A(r) = 
V(r)(^ r ^(r)^'i(r)) is the energy gap [which is assumed to 



be constant (i.e., r-independent) in s-wave superconductors], 
U(r) = -V(r)(*|(r)# t (r)) is the mean (one-electron) po- 
tential, V(r) is a potential characterizing electron-electron in- 
teractions, A(r) is a vector potential, fj, is the chemical poten- 
tial, and a £ {t, I}- Here, (■ ■ • ) implies a combined quantum 
mechanical and statistical average 111 Oil . Note, however, that 
later (to be specified) in this Report, we shall work only at zero 
temperature, wherein the average reduces to (t/>bcs | • ■ • | V'bcs) 
(IV'bcs) 1S the Bardeen-Cooper-Schrieffer state). 

This Hamiltonian can be exactly diagonalized by using the 
following Bogoliubov transformation: 



*t( r ) = E„ u «(r)7» t 
*l( r ) = E„««(r)7ni 



(2) 
(3) 

as H = J2n<r s nlLlna, in which 7^ ( 7nCT ) are 
some bosonic creation (annihilation) operators satisfying 

{lnv,ln><T>} = {lL,ll>a'} = and ilna (?) ,ll> a' )} = 

5aa'5nn'5{Y — r'). Here, £„ is the corresponding energy spec- 
trum, u n (r) and v n (r) are the Bogoliubov coefficients, which 
satisfy 

e n u n {v) = [fle(r) + f7(r)K(r) + A(r)u„(r), (4) 
e n vn(r) = -[H*(t) + £/(r)K(r) + A*(r)u n (r). (5) 

The two-electron space-spin density matrix (TESSDM) — 
up to a normalization factor — in the second-quantized repre- 
sentation is 1151 1 1 ill 



la 1 a 2 ,(j' 1 <j' 2 ( r li r 2! r i: ^2) ~ 



(1-2, 



(6) 



On the other hand, the definition of the two-particle Green 
function is JKH 

G CTl a 2 , (T > iCT > 2 {ri,ti,r2,t2;r' 11 t' 11 r' 2 ,t' 2 ) = (7) 
V[* ffl (ri, ti)* CT2 (r 2 , * 2 )*t ( r ' 2 ,t' 2 )*l, (r'^t',)} 



where T denotes time-ordering, and all field operators have 
been written in the Heisenberg picture. Thus the TESSDM is 
related to the two-particle Green function as follows lIsT Holl : 
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ga 1 a 2 ,a 1 ,a 2 ,(r 1 ,r 2 ;r' 1 ,r' 2 ) = -(1/2) 
x Go-jo-^ojo-/ (ri, ti, r 2 ,t2; r^, tf, r' 2 ,t 2 ), 



(8) 
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where tf = tj + S (with 5 —> + ). Furthermore, the two- 
particle Green function for the superconducting state can be 
written in terms of single-particle Green functions Is UIoll 

G <7l a 2 .cr' 1 cr' 2 (1) 1 , 2') = G CTlCr / (1, 1 )G C r 2Cr ^ (2, 2') 

-G CT1 ^ (1, 2')G CT2ff ; (2, 1') - F airT2 (1, 2)fJ, ct , (1', 2'),(9) 

where the single-particle Green function G(i,j) and anoma- 
lous Green function F(i, j) [here, i = (r,, for the case of 
s-wave superconductors are, respectively, given by ifToll 

G airT[ (n, t x , ri, t[) = -i (T[* ai (ruh)^, (ri, 

= 5 fflCT; G(r 1) t 1 ;r / 1 ,t' 1 ), (10) 

F^ 1CT2 (r l5 ti, r 2j i a ) = -i {T[^ ai (r x , ti)< 2 (r a , i a )]) 

= S 010J ,Ft(r 1 ,ti,r 2 ,t 2 ), (11) 

with 5cr 1(T2 = iuy [recall <r y — (9 ~*J is a Pauli matrix]. 
If we define time-dependent Bogoliubov coefficients u n (r, t) 
and u n (r, t) as in the following equations [everything here, 
except iJ e , has (r, t) dependence, and henceforth we set h = 
c=l]: 



id t u„ = [H e + U]u n + Av n , 
idtvn = -[H* + U]v n + A*u v 



(12) 
(13) 



one can rewrite the Green functions G-|~|-(ri, ii, r^, i^) = 

G(r 1 ,t 1 ,r' 1 ,t[) andF^(r 1 ,t 1 ,r 2 ,t 2 ) = F^(ri, h, r 2 , t<A at 
zero temperature (which we shall assume hereafter), as lUOIl 

G tt (ri,ti,r 2 ,t 2 ) = i^u„(ri,ti)u*(r 2 ,t 2 )f(ii - ^) 

n 

-<(ri,t 1 )«„(r 2 ,t2)fl(*2-*i),(14) 
F^(ri,ti,r 2 ,i 2 ) = i u* (r 2 , t 2 )w„(ri, t\)9(t\ - t 2 ) 

n 

+ «n(r 2 ,<aX(ri,*i)d(t2-*i),(15) 

where 0(f) is the unit step function. 

By assuming that the Hamiltonian and the superconduct- 
ing gap are time-independent, the Bogoliubov coefficients 
can be written as u n (r,t) — w n (r)e _le "* and v n (r,t) = 
w n (r)e~ le "*. We can consider a weak magnetic field as a per- 
turbation, which modifies the Bogoliubov coefficients and the 
energy gap (up to first order) as 

X(r) =X^(r)+X (1 \r)+0(xW), (16) 

in which X 6 {A, u n , v n }, and superscript "(0)" denotes the 
solutions to the unperturbed Bogoliubov equations [Eqs. 
and ©I (i.e., for A = 0). Using the London gauge and con- 
sidering A not having any component perpendicular to the 
surface of the superconducting system, one can ignore A' 1 ' 
ifToh - We can expand u„\r) in terms of a complete set of 
normal-state eigenfunctions, and from Eqs. (|4j, (HJ, and ([Tol l 
[or Eqs. (TTZt and (fT3ll. in which case the linearized Bogoli- 



ubov equations read 



V 2 

£ n + \-(i-U(r) 

Am 



uW(r)-A(°)(r) U W(r) = 



le 
2m~ 



[A(r)-V + V- A(r)]u<, 0) (r) 



(17) 



V 2 

- 2 M + U ( r ) 



U (i)( r )-A(°)*(r)n«(r) = 



le 
2m 



[A(r)-V + V-A(r)]u(°!(r) 



(18) 



We now insert u { „\r) = J2 m a nm <t l rn(r), v£\r) = 

Em & ™™^™( r )' u£\r) = u(£ n )<f>n(T) and ui 0) (r) = 
v(£ n )<p n (r), in which </> n (r) denotes normal-state eigen- 
functions (which can always be taken to be plane waves 
»(r) = e <k "- r / y^), «(£„) = V(l+We n )/2, u(gn) = 
v/(l - ^n/e„)/2, 6> = n 2 /(2m) - fi, e n = ^ + A 2 
(the excitation energy of a quasi-particle with wave vector n), 
and A„ = - J2i V nm u(£ m )v(£ m ) [with the Fourier trans- 
form Km ~ J" V(r)e- i ( n - m )- r d 3 r] Q3. Hence, 

A/f r n 



£2 _£2 
Srt Sri 



, (19) 



(and similarly for b nm by replacing u — >• w and £ m — > — £ m ) 
where 

M„ m = ^ /" [<0«k„ + 0„0;„k m ] • A d 3 r. (20) 

Therefore by calculating the coefficients a nm and b nm , we 
will have m1 (r) and t4* ( r )- From Eq. ( fT5l ), the linearized 
single-particle Green function can be rewritten in terms of 
u n (r) and v n (r), up to first order of perturbation and when 
t 2 —>fi, as follows: 

G(ri,r 2 ) = G(ri,ti,r 2 ,t 2 )| t2 ^ t + = i v* (ri)«„(r 2 ) 

n 

wGWfri.raJ+GWCrx.ra), (21) 

where 

G(°)(r 1; r 2 ) = *( ri )^°)(r 2 ), (22) 

G«( ri ,r 2 ) =^ U ( 1 )*(r 1 >f(r 2 ) +^°)*(r 1 )^ 1) (r 2 ). 

(23) 

More specifically, (the Green function in the absence of 
magnetic field yD) and (the first-order correction) are 
given as 

G( )(r 1; r 2 ) = G( )(r) 



27r 2 r J 



v (Cfc) sin(fcr)fc dfc 



G(°'(0) 



sin(fc^r) 
2k F r 



G (1) (ri,r 2 ) =^&„ m w*(e„)^(r 1 )«^ m (r 2 ) 

+ & nm«(C«) < / , n(r 2 )C( r l)> 



(24) 



,(25) 



(26) 
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FIG. 1. (Color online) Quantum and classical correlations (red, solid 
line; purple, dashed line, respectively) as functions of p of the Cooper 
pair density matrix. Here, entanglement is calculated through con- 
currence (green, solid line), and quantum correlation is calculated 
through quantum discord (which already excludes classical correla- 
tions from quantum mutual information) — Eq. ( 142b . 



in which r = ri — r 2 is the relative distance of electrons 
in a Cooper pair, 1\ (x) is the order 1 modified Bessel func- 
tion, and Li(x) is the order 1 modified Struve function lfl2tl . 
At the continuum limit (k„ — > k, k m — > k', ^2 —> 

[l/(27r) 3 ] J d 3 k), one can obtain 



G«(n,r 2 ) 



in 



Ykk' Mkk< 



(27) 



2(2tt) 6 

X ^(k'^-k-n) +xe i(k-r 2 -k'- ri )] d 3 kd 3 k ' ; 

2 _ C ' 2 \ 



where Y kk , = [(e -${e- £) + A 2 ] /[e(£ 2 - O], M kk , is 
the continuum version of Eq. (l20l i. and \ is +1 (—1) if M kk i 
is real (imaginary). 



Additionally, we have 

F\v 1 ,y 2 ) = F\ri,t 1 ,r 2 ,t 2 )\ t2 ^ t + = i ^ u n (r 2 )<(ri) 

n 

wF(°»(r 1) r 2 ) + ^)t(r 1)rj ) ) (28) 

where 

F(°)t( ri ,r 2 ) =i^2v^(r 2 )u^*(v 1 ) (29) 

n 

F( 1 )t(r 1 ,r 2 )= l ^ u W(r 2 ) U (°)*(r 1 )+^)(r 2 ) U W*(r 1 ), 

n 

(30) 

or equivalently 

Ft(°)(r 1 ,r 2 ) = i^v(^)^(r 2 K(^)^(r 1 ) 

n 

= p(0)*^^ _ * f sin(fcr)fc 



F w *(r) = 
iAiV(O) 



sin(fci?r) 



dfe 



W-fY 



(31) 



F«t( ri) r 2 ) = i^6„ mU *(e„)C(ri)0 m (r 2 ) 



+ a™m«(Cn)0n( r 2)</>m(ri), (32) 

where iV(0) = (1/2) [fc 2 (dfc/de fc )] = mfc F /(27r 2 ) is 

the density of states for one spin projection at the Fermi sur- 
face |Ht], and Ko(y) is the order modified Bessel function 
fl2ll . and in Eq. 07\ and here A is A^ (constant for s-wave 
superconductors). At the continuum limit, one can also obtain 



F (1)t (ri,r 2 ) = 
1 



iflA 



2(2tt) 6 



d 3 k d 3 k'^^{e- 4 ( k - ri - k '- r2 > 



i(k-T2— k -ri) 



}■ 

(33) 



The TESSD, for the case (ri , r 2 ) = (r[ , r 2 ) and up to first 
order approximation, reads J5[ 



£W 2 , CT ^(ri,r 2 ) = -(1/2) S aia[ S a ^ 2 {G^ 2 (0) + G ( °Ho)G^ (r 2 ,r 2 ) + G (0) (0)G (1) (ri, n)} 
-6^4 <Wi {G(°) 2 (r) + G<°> (r)G« (r 2 , n) + G<°> (r)G« (n , r 2 )} 
-^ lCT2 ^ ; ^{|F (0) M| 2 +^ (0) M^ (1) *(ri,r 2 ) + F(°)*(r)^ 1 )(r 1 ,r 2 )} 



(34) 



Using n(r) = — iTr[G(r, t, r, t + )], we can write 
G^^(ri,ri) = — m(ri)/2 (similarly for n(r 2 )) and 
G(°)(0) = -in/2, where n is the electron density. We thus 
obtain 



/Pi 



Tr[g] 2{ Pl +p 2 ) 



P2 P3 
P3 Pi 



(35) 



Pi 



where 



pi =d + q+(l-g 2 )n 2 /8, 
p 2 = d + w+(l + f)n 2 /8, 
p 3 = q-w-(g 2 + f 2 )n 2 /i 



(36) 
(37) 
(38) 



with/(r) = F(°)(r)/G(°)(0), g{r) = G<®(r)/G<®(p) 
2iG(°)(r)/niU(r l7 r 2 ) = [n( ri ) + n(r 2 )]n/8, g(r Xj r 2 ) 
-i[GW(r 2 ,n) + G( 1 )(r 1 ,r 2 )] 9 (r)n/4, and W (r 1; r 2 ) 
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Re[F(°)(r)f (1) *(ri,r 2 )] 
g is in fact a Werner state 



13|j. It is straightforward to see that 
IB 



1-p 



l+p\if>-){il> |, 



(39) 



where 1 is the 4 x 4 identity matrix, \ip ) = (| f-J-) - | 4-t 
))/v2 is the spin singlet, and 



w - 9 + (/ 2 + S 2 )n 2 



2d + w + q + (f 2 - g 2 + 2)n 2 / 



(40) 



It is interesting to see how the state (and the correlations 
thereof) vary due to the application of the magnetic field (cf. 
Ref. H]). Neglecting f(r) 153 and terms higher than 0(G (1) ), 
the perturbed p [denoted below as p^> ] has the following form 
in terms of the unperturbed one which is for the case without 
magnetic field [denoted below as p(°']: 



) g 2 (r)d(ri,r 2 ) + gQi,r 2 ) 
[2-3 2 (r)] 2 



(41) 



As a result, it can be argued that for some values of the 



magnetic field, p 



(i) 



P 



(0) 



This implies that the addition of 



a nonvanishing, weak magnetic field sometimes does not alter 
the related density matrix (and hence, as we see below, its 
quantum correlations), while in general it affects the density 
matrix as in Eq. d4"TT >. 

Having obtained the form of the density matrix, it is usually 
of interest to know how it is correlated quantum mechanically. 
To that aim, entanglement is usually considered as a measure 
of quantum correlation. Using the properties of a Werner state 
flql . it is straightforward to calculate entanglement through, 
for example, the Peres-Horodecki separability criterion H13fl - 
Specifically j_a 2x2 Werner state is entangled if and only if 

p > 1/3 sum]. 

Although non-separability (i.e., entanglement) is consid- 
ered to be a purely quantum mechanical effect, it has been 
known that it is not responsible for all quantum correlations. 



To quantify non-classicality of correlations, "quantum dis- 
cord" has been introduced as a finer measure. It is defined 
to be the difference between two classically identical expres- 
sions for the mutual information 01 ■ While no separable state 
can have quantum entanglement, some of them may exhibit 
nonvanishing quantum discord. For Werner states of the form 
d39l ). it has been shown that any p > implies a nonvanish- 
ing discord |7lll6t]. More explicitly, quantum discord can be 
obtained with the following relation flol : 



Q(q) 



1 r 

4 



(1 - p) log 2 (l - p) + (1 - 3p) log 2 (l + 3p) 
2(l+p)log 2 (l+p)l, (42) 



which can be expanded based on Eq. (l4"TT i to see the effect of 
the magnetic field. Figure[T]depicts quantum discord and con- 
currence ifnll (as a specific measure of entanglement, equal to 
max{0, (3p — l)/2}) for the spin state of a Cooper pair. 

In summary, we considered an s-wave superconductor sub- 
jected to a weak magnetic field, and obtained spin-space den- 
sity matrix of electrons of a Cooper pair. We employed the 
relation of the density matrix with the Green functions, and 
applied perturbation theory to solve the corresponding Bogoli- 
ubov equations. This density matrix was found to be in the 
form of a Werner state, hence it is not necessarily maximally 
entangled. This can be attributed to the fact that correlations 
between electrons in the ground state of a superconducting 
state are not localized in each Cooper pair |5|]. Nonetheless, 
a nonvanishing quantum correlation may have some relevant 
implications on the overall quantum correlation in the system. 
We thus calculated quantum discord — as a measure of quan- 
tum correlation — for the obtained density matrix, and showed 
how a weak magnetic field affects this correlation. Besides 
basic importance of such calculations, we hope that our study 
can shed some light on potential applications of supercon- 
ductors in, for example, quantum information science through 
nanoscale systems itlal . 
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